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Abstract 



The quartic Henon-Heiles Hamiltonian H = (Pf + P£ )/2 + (OiQf + VL 2 Qi)/2 + CQf + BQfQi + AQ\ + 
(l/2)(a/Qf + P/Q2) — 7Q1 passes the Painleve test for only four sets of values of the constants. Only 
one of these, identical to the traveling wave reduction of the Manakov system, has been explicitly 
integrated (Wojciechowski, 1985), while the three others are not yet integrated in the generic case 
(a, /3, 7) 7^ (0, 0, 0). We integrate them by building a birational transformation to two fourth order first 
degree equations in the classification (Cosgrove, 2000) of such polynomial equations which possess the 
Painleve property. This transformation involves the stationary reduction of various partial differential 
equations (PDEs). The result is the same as for the three cubic Henon-Heiles Hamiltonians, namely, 
in all four quartic general solution which is meromorphic and hyperelliptic with genus two. 

As a consequence, no additional autonomous term can be added to either the cubic or the quartic 
Hamiltonians without destroying the Painleve integrability (completeness property). 
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1 Introduction 



The considered Hamiltonian originates from celestial mechanics, as a system describing the motion of a 
star in the axisymmetric potential of the galaxy. Denoting q\ the radius and qi the altitude, this "Henon- 
Heiles Hamiltonian" (HH) JS| is the sum of a kinetic energy and a potential energy, in which the potential 
is a cubic polynomial in the position variables q\ , qi , 

H = \{p\+pl + ql + ql) + qiql-\ql, (1) 

it is nonintegrable and displays a strange attractor. However, if one changes the numerical coefficients in 
the potential, the system may become integrable, and this question (to find all the integrable cases and to 
integrate them) has attracted a lot of activity in the last three decennia. 

A prerequisite is to define the word integrability ; and in section [3 we briefly recall its three main accep- 
tations in the context of Hamiltonian systems. 

In section [21 we recall all the cases (three "cubic" plus four "quartic" ) for which the most general 
two-degree of freedom classical time-independent Hamiltonian may have a single valued general solution. 
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Then, discarding the integrated cases (see f° r a review of the current state of this problem), we focus 
on the three cases (all "quartic" ) for which the general solution is still missing, with the aim of finding this 
general solution. 

In section we build an equivalent fourth order ordinary differential equation (ODE) for qi(t), in the 
hope of finding it listed in one of the classical tables of explicitly integrated ODEs. This hope is deceived 
because these tables are not yet finished. 

This is why, in the last two sections, we adopt a different strategy. In front of the difficulty to perform the 
separation of variables in the sense of Arnol'd and Liouville, we establish a birational transformation between 
the two second order Hamilton equations and a fourth order ODE listed in a classical table established by 
Cosgrove [TT], whose general solution is single valued. 



2 Integrability for Hamiltonian systems 

Given a Hamiltonian system with a finite number N of degrees of freedom, three main definitions of 
integrability are known, 

1. the one in the sense of Liouville, that is the existence of N independent invariants Kj whose pairwise 
Poisson brackets vanish, {Kj, K{\ = 0, 

2. the one in the sense of Arnol'd-Liouville 2 , chap. 9] , which is to find explicitly some canonical variables 
Sj, Tj,j = 1, N which "separate" the Hamilton- J acobi equation for the action S, which for two degrees 
of freedom writes as, 

dS dS 

H(qx,q 2 ,pi,p 2 ) - E = 0, p x = — , p 2 = ^— , (2) 

oqi oq 2 

3. the one in the sense of Painleve |S] i.e. the representation of the general solution qj(t) by an explicit, 
closed form, single valued expression of the time t. 



3 The seven Henon-Heiles Hamiltonians 

Given the most general two-degree of freedom classical time-independent Hamiltonian 

H = \{j>l+pl) + V{ qi ,q 2 )=E, (3) 

the requirement that the system made of the two Hamilton equations passes the Painleve test [Hj (for at 
least some integer powers q™ 1 , q% 2 ) selects seven and only seven potentials V depending on a finite number 
of constants, namely 

1. three "cubic" potentials (HH3 case) |7lH31l9^. 

H = ^{pI+P 2 2 + wig? + Lo 2 ql) + a qi ql - -f3ql + \iq 2 2 , a^Q (4) 

in which the constants a, (3, oj\, uj 2 , 7 can only take three sets of values, 

(SK): / 8/a=-l,«i = tJ a , (5) 
(KdV5) : P/ a = -6, (6) 
(KK) : (3/a = -lG, w 1 = l6w 2 . (7) 
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2. four "quartic" potentials (HH4 case) |21IT7| . 

1 



H = ^(Pi + Pi + fiiQi + fiaQl) + CQf + BQlQl + AQ 4 2 



in which the constants A, B, C, a, f3, 7, ili,fl 2 can only take the four values (the notation A : B : C = 
p : q : r stands for A/p^B/q = C/r — arbitrary), 

A : B : C= 1 : 2 : 1, 7 = 0, 

A:B: C= 1:6:1, 7 = 0, Q X = fi 2 , , , 

A : B : C = 1 : 6 : 8, a = 0, O x = 4ft 2 , l j 

A : B : C = 1 : 12 : 16, 7 = 0, fii = 4f2 2 - 

All seven cases are integrable in the sense of Liouville, with a second constant of the motion K [121 ^5] 
[201 El E] either quadratic or quartic in the momenta Pi,p%. 

In the sense of Arnol'd-Liouville, the separation of variables has been performed [T21 13T1 051 051 l6l OS ] . 

except in three cases, 

1. HH4 1:6:1 a ^ (3, 

2. HH4 1:6:8 /3 7 ^ 0, 

3. HH4 1:12:16 af3 ^ 0. 

What is remarkable is the fact that, in all cases when the separation of variables is achieved, the equations 
of Hamilton have the Painleve property, the general solution being a hyperclliptic function of genus two. 
The purpose of this work is to prove equally the Painleve property in the three remaining cases where the 
separation of variables is not yet performed. 

4 Equivalent fourth order ODEs 

In the cubic case, the two Hamilton equations 

q';+uj iqi -Pql + aq* = 0, (10) 
q 2 + u 2 q 2 + 2aq x q 2 - jq 2 3 = 0, (11) 

together with the Hamiltonian are equivalent ^5] to a single fourth order ODE for qi(t), 

20 

ql" + (8a - 2/3) qi q1 - 2(a + (3)q' 2 - -a^q 3 

+(u>! + 4u 2 )q" + (6auj 1 - 4:j3uj 2 )q 2 + 4cjicj 2 «i + 4aB = 0, (12) 

independent of the coefficient 7 of the nonpolynomial term q^ 2 and depending on the constant value E of 
the Hamiltonian H. In the three HH3 cases 10-0, this ODE belongs to a list JI] ("classification") of 
equations enjoying the Painleve property, whose general solution is hyperelliptic with genus two. 

In the quartic case, the similar fourth order equation is built by eliminating Q 2 and Q'{' 2 between the 
two Hamilton equations, 



Q'l + fiiQi + ACQl + 2BQ 1 Q 2 - aQf + 7 = 0, (13) 
Q 2 + fl 2 Q 2 + 4AQ 3 + 2BQ 2 Q 2 - f3Q^ 3 = 0, (14) 
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and the Hamiltonian 



which results in 
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(15) 



This ODE depends on E but not on (3 and, as opposed to the cubic case, it does not belong to a classified 
set of equations, because Q'{" is not polynomial in Q\. 

In the three remaining cases, since one is yet unable either to perform the separation of variables or to 
establish a direct link to a classified ODE, let us build an indirect link to such a classified ODE. This link, 
which involves soliton equations, is the following. 

For each of the seven cases, the two Hamilton equations are equivalent [151 1161 13"] to the traveling wave 
reduction of a soliton system made either of a single PDE (HH3) or of two coupled PDEs (HH4), most 
of them appearing in lists established from group theory |13| . Among the various soliton equations which 
are equivalent to them via a Backlund transformation, some of them admit a traveling wave reduction to a 
classified ODE. This property defines a path [22U3U) which starts from one of the three remaining HH4 cases, 
goes up to a soliton system of two coupled 1+1-dimensional PDEs admitting a reduction to the considered 
case, then goes to another 1+1-dim PDE system equivalent under a Backlund transformation, finally goes 
down by reduction to an already integrated ODE or system of ODEs. 



5 General solution of the quartic 1:6:1 and 1:6:8 cases 



Let us denote the two constants of the motion of the 1:6:1 and 1:6:8 cases as, 

1 



H 



1 : 6 : M 



and 



1 : 6 



K 



K = 
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2 -**> 
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2 \Q\ Q2 
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(Qt + 6QiQi + Q\ 



Ql 
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E, 
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-.2^.2 



TiQV 



H=-(pt+ Pi) + ~(4<rf + <g) - —(8qt + 6qfqi + q^ 



+ 7 ( -27172 - 4q 2 p!P2 + -q\q\ + q\ql + ^q\p\ - 4^(?i<72 + 
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(16) 



(17) 
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There is a canonical transformation [5] between the 1:6:1 and 1:6:8 cases, mapping the constants as follows, 

Ei-.e-.s = Ei :6 a, ^1:6:8 = ^1:6:1, W = O, 7 = M - — , /3 = -(«! - K 2 ) 2 , (18) 

therefore one only needs to integrate either case. 

The path to an integrated ODE comprises the following three segments. 

The coordinate qi(t) of the 1:6:8 case can be identified 0J|3j to the component F of the traveling wave 
reduction f(x,r) — F(x — cr), g(x,r) = G(x — cr) of a soliton system of two coupled KdV-like equations 
(c-KdV system) denoted c-KdVi glEj 

fr + + \ff x - \f + = 0, 

-2ffr + g X xx + 6gg x + 3fg xx + 6gf xx + 9f x g x - 3f g 

(19) 

33 22 

^fxxxx ^ffxxx ~l~ Qfxfxx 3/ f xx 3/ f x = 0, 

with the identification 

qi = F, q\ = -2 (F' + F 2 + 2G- 2w) , , * 

c = -w, K± = 7, K 2 = E, [ ' 

in which K\ and -ft^ are two constants of integration. 

There exists a Backlund transformation between this soliton system and another one of the c-KdV type, 
denoted bi-SH system [T3"| . 

-2u r + (u xx + u 2 + 6v) x = 0, , , 

+ v xxx + uv x = 0. 

This BT is defined by the Miura transformation [221 

«=§ (2.g -/.-/ 2 ), 

3 (22) 
r - (2/,xx + 4// xx + 8 5 /, + 4/ 5;c + 3/ 2 - 2f 2 f x /' • 4g/ 2 ) . 

Finally, the traveling wave reduction u[x, r) = {/(a; — cr), u(:c, t) = V^o; — cr) can be identified | 30| to 
the autonomous F-VI equation (a-F-VI) in the classification of Cosgrove jll| . 

2 

a-F-VI : y"" = 18t/y" + 9t/ 2 - 24y 3 + a vl y 2 + ^-y + Kyl t + fa, Kyl = 0, (23) 

an ODE whose general solution is meromorphic, expressed with genus two hyperclliptic functions [111 
Eq. (7.26)]. The identification is 

u = ~ 8 (» + TS. 

V = y"~6y 2 + ^cy+^-c 2 -^, (24) 

512 

a V i = -4c, fa = K B - 2cK a + ^c 3 , 

in which Ka,Kb are two constants of integration. 

In order to perform the integration of both the 1:6:1 and the 1:6:8 cases, it is sufficient to express (F, Q) 
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rationally in terms of (U, V, U', V). The result is 



F 



W 



2W 24W 



-W 2 - 2{U - 3c) (12V +(U + 3c) 2 ) 
+ 36K B - 54K 2 ] , 



G = f + ^ 2V + 3K2) ( 2V " + KlU '~ 3K i) 

- 2(U - 3c) (2K 1 V + Kf(U + 3c))] , 

2 



(25) 



W 



V + -K, ) +^ 1 2 (C/-3c), 



I K A = K 2 . 

Making the product of the successive transformations (|20|l . I|25|l. I)24[l. one obtains a meromorphic general 
solution for Q 2 , Q 2 ., <7ii 92 > 



9i 



9l 



W 7 
2W + W 



9j-3[y+^)(h + E)-^ 
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-16 ^ - -a; 
1 

W 



12 (V + ^ V - 48y 3 - 16wy 2 + ( 24£ + i^cj 2 ) y + 
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1280 
243 L 



- ^loE + 24 7 (y - fc' - 144 7 2 (y - ~w 

W=(h + E) 2 - 9 7 2 (y - ^uA , 

3 3 512 

a V i = 4w, /?vi = 7 7 2 + 2^ - — /3 - — w 3 , 

3 l 4 3 6 24 ? 9 

^ = 32* - 2*> * 2 < VI = 32** " 3* + 64^< 

#i = 7, #2 = -E, = E, K B = -— /3 + - 7 2 . 

16 4 

in which /i and j are convenient auxiliary variables |111 Eqs. (7.4)-(7.5)], 

Q{sx,s 2 ) + y/Q(s 1 )Q{s 2 ) 5 



36 



3^ Sl s 2 + Cvi + V( s i - Cvi)(s* - Cvi) F V i 
ft — --tiyi — . 



1 



j = _(2/i + F V i)<iy + 



Sl + S2 

«vi £vi 



9 4(si + s 2 ) 



In the above, the variables s\, s 2 are defined by the hyperelliptic system 11 



(si - S2K1 = VP(si), («2 - si)s 2 = a/P(s 2 ), 
P(s) = (s 2 -Cvi)Q(s), 

Q( S , t) = (.s 2 - C V i)(t 2 - Cvi) - + t 2 - 2C V i) + ^f-(s + t) + F Vh 

. Q(s) = Q{s, s). 



(26) 



(27) 



(28) 



The expressions (|27|l cannot be written as rational functions of s±, s 2 , s[, s' 2 and are nevertheless meromorphic 

[HEP 

The coefficients (a, C, E, F)vi of the hyperelliptic curve depend algebraically on the parameters of the 
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Hamiltonians P,j,u,E,K HU Eqs. (7.9)-(7.12)] 



A V i = 4uj, 

= -yc^vi -2E)-0 + 4 7 2 , 

(Fvi - 2E) 2 (F VI + AE) + — (Fyi - 



(29) 



2E) 



27 



/3 7 2 = 0, 



and this algebraic dependence could explain the difficulty to separate the variables in the Hamilton- Jacobi 
equation. Note that, in the particular case (3-f = 0, i.e. n\ = «|, these coefficients become rational, see |29| . 
Remark. The F-VI ODE can be written in Hamiltonian form, 



H 



Qj 3 
6 2 



13 11 

34T a 9 



2 9 3 3 



VI 



(30) 



^ Qi = -6 (y - ^) , Q 2 = -62/', P x = 6?/'" - 108w', P 2 = -6y" 



In the autonomous case kvi = 0, the Hamiltonian H is a first integral (equal to 36X1^1), and the other 
constant of the motion is cubic in the momenta. However, because of the nonlinear link between Kiyi 
and the two first integrals of the 1:6:8 case, see (|24|l . there exists no canonical transformation between the 
variables (qj,Pj) of 1:6:8 and the above canonical variables of a-F-VI. 

6 General solution of the quartic 1:12:16 case 

Let us denote the two constants in involution as, 



h = \{P 2 , + Pi) + %Qi + Qi) - ^(i6Qt + UQlQl + Qi) 

1 f k\ 4k| 

T2 J 



E, 



2 V Qi Qn 

1 : 12 : 16 < i / K a 

K = — \ 8(Q 2 P - QiP 2 )P 2 - QiQ 2 - 2Q\Ql + 2ilQ 1 Q 2 2 + 32Q 1? | 

16 v ~ Qi 



Ql 



(31) 



Similarly to the 1:6:1-1:6:8 couple, there exists a canonical transformation between the 1:12:16 Hamil- 
tonian and another Hamiltonian [3J 0] , which is however not the sum of a kinetic energy and a potential 
energy, which we denote similarly as 5:9:4, 



5 : 9 : 4 I 



H 



Pi - 29i<72 



- \{K + Hol + 5<Z 2 4 ) + |(«? + ql) - KQ! + ^ = E, 



1 2 
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9 2 
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(32) 
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The path to an integrated ODE is also quite similar and is made of the following three segments pi l27l!22| . 

Firstly, the coordinate qi(t) of ilML'l is identified 0J|3] to the component F of the traveling wave reduction 
f(x,r) = F(x — cr),g(x,T) = G(x — cr) of a soliton system of two coupled KdV-like equations denoted 
c-KdVa(/,<?) HE], 



F, q 2 2 = -(F'-2F 2 -G 



(33) 



Secondly, there exists a Backlund transformation between this soliton system and another one of the 
c-KdV type, denoted bi-SK system (u, v) |23|. transformation defined by the Miura map 



i - ^ (3/x - f + 2g) , 

7 = To ^ xxx + 9xx + ^ x9 _ ^ 9x ~ ^ xx + g2 ^ ' 

Finally, the traveling wave reduction u(x, r) = U(x — cr), r) = — cr) is identified j23], 



(34) 



U ■■ 



V = -6y" + I8y 2 



—uj z - -E, 
10 5 ' 



(35) 



to the F-IV equation (or to the F-III as well) in the classification of Cosgrove ^I] j 



y"" = 30yy" - 60y 3 + a w y + fry, 

1 / / / r> 
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IV 



S 2 )s'i = \/P{8i), (82 ~ Si) S ' 2 = y/P(s 2 ), 

Ay 

3 v ~ ■ - ■ - 3 
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(36) 



3B 



G4 



in which (KijVi ^2,iv) denote two polynomial first integrals of F-IV. The general solution of this ODE is 
meromorphic, expressed with genus two hyperelliptic functions 

In order to perform the integration of both Hamiltonians (|31|) and l|32|) , it is sufficient to express (F, G) 
rationally in terms of (U, V, U', V). The result is 



w 

F = +K X a X 2 , 

2W ' 

G = -F 2 - X X X 2 + Ki 
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X 1 = V + 2U' - ZuU + —uj z E, 
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25 2 V 1 2J 2 h ' 
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25 



27 
T25 C 



From the point of view of the separation of variables, one should first exhibit a Hamiltonian represen- 
tation of F-IV. One such structure is that of the cubic SK case. However, since the constant value of the 
Hamiltonian of the cubic SK case, when expressed only in terms of the parameters (E, K, oj, n\ 1 k 2 ) of the 
1:12:16, is not an affine function of E, there exists no canonical transformation between the cubic SK case 
and the 1:12:16 case. 
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7 Conclusion, remaining work 



The explicit integration of all the seven cases is now achieved in the Painleve sense (finding a closed form 
single valued expression for the general solution), and the common features are the following. 

1. In all cases, the general solution is hyperelliptic with genus two, and therefore meromorphic. 

2. Each case is birationally equivalent to a fourth order ODE which is complete in the Painleve sense, 
i.e. which accepts no additional term, under penalty of losing its Painleve property. Consequently, for 
each of the seven Hamiltonians, it is impossible to add any term to the Hamiltonian without destroying 
the Painleve property, and the seven Henon-Heiles Hamiltonians are complete. 

About the integration in the Arnol'd-Liouville sense (finding the separating variables of the Hamilton- 
Jacobi equation), two problems remain open. 

1. In the 1:6:1-1:6:8 case, the hyperelliptic curve y 2 = P(s) of F-VI (see 1)28(1 1 reduces in the separated 
cases /?7 = to the hyperelliptic curve of the separating variables. Therefore, F-VI is the good ODE 
to consider, and the only missing item is to find a Hamiltonian structure of F-VI, necessarily distinct 
from (|3(jp. admitting a canonical transformation to 1:6:1-1:6:8. 

2. In the 1:12:16-5:9:4 case, the hyperelliptic curve y 2 = P(s) of F-IV (see (jHEJ)) does not reduce in the 
separated cases K1K2 = to the hyperelliptic curve of the separating variables, which is 

Kin 2 = : P(s) = s 6 - lus 3 + 2Es 2 + — s + k\ + k\ = 0. (38) 

Therefore, F-IV (as well as its birationally equivalent ODE F-III) is not the good ODE to consider, 
and it should be quite instructive to integrate the fourth order equivalent ODE (|15fl in that case. 
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